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Abstract 

We employ a molecular model to study a ternary mixture of saturated lipid, with tails of sixteen 
carbons, a mono unsaturated lipid with tails of eighteen carbons, and cholesterol. The model, 
solved within mean-field theory, produces several forms of phase diagrams depending upon the 
relative strengths of interactions, but only one that shows the coexistence of two liquid phases 
observed in experiment. The lipids in the phase rich in cholesterol are more ordered than those in 
the other. The binary cholesterol, saturated lipid system also exhibits liquid, liquid coexistence. 

PACS numbers: 64.70. Ja,87.16Dg,81.30Dz 
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There has been enormous interest in the hypothesis that the hpids comprising the plasma 
membrane are distributed inhomogeneously, with regions rich in cholesterol and saturated 
lipids floating, like rafts, in a sea of unsaturated lipids. Such rafts have been implicated in 
many biological processes including endocytosis, transcription and transduction processes, 
and viral infection. The extent of the interest is reflected in the number of recent reviews 
^lO^y^y^Ol- It has spurred numerous in vitro investigation of ternary systems of cholesterol 
and a high-melting point (hmp) lipid, either sphingomyelin or a saturated glycerolipid, and 
a glycerolipid whose melting point is lower, either because it is unsaturated 0, 0,0, m 
because, .s „ major question the ternary 

studies could address is whether the aggregation phenomenon indicated in vivo might simply 
be attributed to liquid-liquid phase separation. If so, it could be observed in vitro. The 
assumption was that one fluid, the liquid disordered (Id) phase, would be rich in the low- 
melting-point (Imp) lipid, while the other, the liquid ordered (lo) phase, would be rich in 
cholesterol and the hmp lipid. The nomenclature "liquid ordered" I4I signifies that there 
are fewer thermally excited gauche bonds in the hydrocarbon tails of the lipid chains in 
this phase than in the liquid disordered phase. In addition to these two phases, a gel phase 
occurs below the melting, or main chain transition temperature, of the hmp lipid when the 
concentration of this lipid is large. This phase is the most ordered of the three. 

While there is some agreement on broad features in the interior of the ternary diagrams, 
there is also considerable uncertainty. In particular the approach to the cholesterol, hmp lipid 
binary system, and the phase diagram of that system itself remain a matter of controversy. 
The disagreement centers on whether or not there is coexistence of two liquid phases over 
some range of temperature even in this binary system. 

Evidence of such coexistence derives from several studies using different probes: electron 
spin resonance ilSl, freeze fracture electron microscopy llJI, nuclear magnetic resonance 

n n n 

(NMR) 115J, fluorescence resonance energy transfer and other fluorescence studies 
Perhaps the strongest evidence, however, is indirect, coming from calorimetry 

HQ. In 

the cholesterol, dipalmitoylphosphatidylcholine (DPPC) system, these measurements clearly 
show a very large and narrow specific heat signal at a temperature which is rather constant 
over a range of compositions, the standard indicator of a triple point. The signal was so 
interpreted by Vist and Davis [15]. Given the undisputed existence of a gel phase which 
coexists with a liquid phase at low temperature, the observation of a triple point leads to 
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the conclusion of the existence of a gel and two liquid phases at higher temperatures. There 
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is NMR evidence for such a triple point in another cholesterol, hmp lipid system as well 

There are also arguments against the existence of two liquid phases in the cholesterol, 
hmp lipid system. First, the sharp specific heat signals have been interpreted as not arising 
from three-phase coexistence by McElhaney and co-workers [17], and they also dispute the 
conclusion drawn by Vist and Davis from their NMR experiments. Indeed the NMR work 
of Huang et al. |l9| found no such coexistence, nor did those using fluorescence recovery 
after photo bleaching |2o| . Fluorescence microscopy also does not observe phase separation 
in these systems 0, 0|- Even if the behavior of this binary system were clear, its relation 
to that of the canonical 1:1:1 raft composition would not be. The question would remain 
whether the coexistence of liquids in this part of the ternary system was directly related to 
such coexistence in the binary system. 

To clarify the ternary system and its evolution from the binary one, we have carried out a 



study based upon a microscopic model. Previous theoretical investigations of t 



lese mixtures 



have been limited to the binary systems of cholesterol and hmp lipid or hmp 

and Imp lipids j^] • Notable among the formerjs the pioneering study of Ipsen et al. which 



coined the description "liquid-ordered phase" 



121. 



Our ternary mixture consists of cholesterol, a lipid with two saturated tails of sixteen 
carbons, C16:0, and another with two mono-unsaturated tails of eighteen carbons, C18:l. 
There are essentially two adjustable interaction parameters in the model, and their variation 
produces three distinct phase diagrams. One is clearly irrelevant to experiment, one is con- 
sistent with claims of no liquid, liquid coexistence in the binary system, and one is consistent 
with claims that there is. Only the last of these exhibits the liquid, liquid coexistence in the 
heart of the ternary diagram, consistent with experiments. 

In our model, the lipids tails are described using the Flory rotational isomeric states 
representation [2J|. One introduces the local volume fraction per unit length for a lipid tail 
in a particular configuration, (pa{z), 

4><t{z) = ^l'a{k)S[z - Zk\, (1) 
fc=l 

where the z direction is perpendicular to the plane of the bilayer, a = s or u is an index 
indicating whether the lipid is saturated or unsaturated, the z/(A;) = 28A^ = z/q, the volume 
of the /c'th monomer if it is a CH2, or 56A^ if a CH3 group, and = 15, = 17. The local 



orientation of the chain is conveniently specified by the normal to the plane determined by 
the /c'th CH2 group, 

= 1 r, «; = ^--Us - 1. (2) 

A similar local volume fraction is introduced for the cholesterol, with index a = c, ric = 27 
and z/c(A;) ~21.0A^. The orientation of the small acyl chain of the cholesterol is specified in 
the same manner as are the lipid chains, while the orientation of its rigid rings are specified 
by the unit vector, Uc, from the third to the seventeenth carbon in the molecule, using the 
conventional labeling. 

In ea.„e. theo.es. QQ . the effeet of the sho.t.ange .epuls.ve and long-.ange attractive 

interactions between elements was accounted for approximately by replacing them by a 
constraint that the density within the hydrophobic region be constant locally. The free 
energy of the system can then be written in terms of the 0. While accounting well for many 
properties of bilayers, this theory cannot lead to a main chain transition to a gel phase. 



Hence we have added 



23| an additional separable interaction per unit volume which tends 



to align the elements with each other and with the bilayer normal c: 

= -(Ja,a'/'^0)5'(UCT ■ c)5f(uCT/ ■ C), (3) 

(7(u.-c) = (m + l/2)(u.-c)2- (4) 

For large m, (? ~ mexp(— m^^) where Q is the angle between the two unit vectors. Matching 
lipid parameters, we have taken m = 18. Note that this interaction between two vectors 
falls off exponentially if either of them is not well aligned with the bilayer normal. We 
take the strength of the local interactions between bonds in lipid tails to be the same, 
Jss = Juu = Jsu = Jih irrespective of whether the bonds are in a saturated or unsaturated 
chain, and the interaction strength between the cholesterol and any bond in a lipid tail to be 
the same, Jgc = Juc = Jic, irrespective of the type of lipid chain. Thus the model contains 
three interaction strengths: Ju which is set by the main chain transition temperature, Jic, 
and Jcc the strength of the aligning interaction between cholesterols. 

Just as the describe the local density, so we introduce a function which describes the 
local ordering with respect to the bilayer normal 

rto- — 1 

L{^) = J2 ^cT{k)5{z - Zk)g{\ia ■ c). (5) 

k=l 
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The Helmholtz free energy per unit area, fA(T,ps,Pu,Pc), with p^- = N^j/A, cr = s or u, 
the areal density of saturated or unsaturated chains, and pc = N^/A the areal density of 
cholesterol, can now be obtained directly within mean-field theory. One finds 

f^fA =-— /{(/3/2) J--'P- < > < > 

+ Y.P'^[<^-> Ba{z)+ <k>^{z)]}dz 
u 

- V In Q,^ + ^ In — + ^ In — + Pc In — (6) 

a 2 Pq 2 Po Po 

with Po = V/Avq, (3 = l/ksT. Here < Oa{z) >= TrPf^Oa{z), the average in a single chain 
ensemble, with 

Pa = -^exp{-£:^} = -^expi-lSHi 

- - I [^a{zMz) + Uz)Ba{z)]dz}, (7) 
= Trexp{-S„}, (8) 

and Hi contains the intra chain energy arising from the presence of gauche bonds. The 
three unknown fields Bg = Bu = Bi, B^, and 11 are obtained from minimization of the free 
energy with respect to the < ia > and < (p^ > which leads to the three local self consistent 
equations 

I = Y.p,<^^{z)>, (9) 

cr 

Bi = -f3{Jii[ps < L > +pu < L >] + Jicpc <L>}, 

Be = -P{Jlc[ps <L> +Pu < L >] + JccPc < L >}■ 

The first of these is simply the incompressibility constraint. 

The heart of the above method, and its difficulty, is the evaluation of the single molecule 
partition functions Q^j. To evaluate them, we have generated on the order of 10^ configura- 
tions of each molecule. Finally the effect of the water, lipid interface is taken into account 
via a contribution to the free energy per unit area of 70, set equal to the oil, water tension. 
Thus the total free energy per unit area, or surface tension, is ftotiT, ps, pu, Pc) = /a + 7o- 
We look for phases for which this surface tension vanishes. Phase equilibria is determined 
by standard thermodynamic equalities [2^ . 

We now turn to the results of our calculation. The binary phase diagram of unsaturated 



and saturated lipids 



23| depends upon packing constraints and the lipid-lipid interaction 
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FIG. 1: Calculated phase diagrams of the three binary mixtures of cholesterol (c), saturated (s), 
and unsaturated (u) lipids in temperature-composition space for Jii{m + l/2f/kBT* = 1.44 and 
Jic = Jcc = 0.0. These binary diagrams form the sides of the Gibbs prism, a cut through which at 
300K produces the Gibbs triangle shown in Fig. lb. Regions of two-phase coexistence are shaded, 
and some tie lines are shown. 

strength Ju which we set by matching the main chain transition of the saturated hpid to 
T* = 3157^, that of DPPC Below this temperature there is a gel phase, rich in the saturated 
lipid, and a disordered liquid phase rich in the unsaturated lipid. The main chain transi- 
tion temperature of the unsaturated lipid is below 0°C. The nature of the ternary diagram 
depends upon the relative strengths, J/c and Jcc of the lipid-cholesterol and cholesterol- 
cholesterol aligning interactions. We find essentially three classes of diagrams. 

The first occurs when the other two interactions are weak compared to Ju. For the binary 
cholesterol, saturated lipid system, this results only in gel and Id phases, and therefore no 
liquid, liquid coexistence, consistent with one side of the controversy noted earlier. The 
three binary phase diagrams as a function of temperature and composition are shown in 
Fig. 1(a). A ternary diagram at a temperature below the main chain temperature is shown 
in Fig. 1(b). There is no region of liquid-liquid coexistence. As such regions can be seen 
directly in the experimental ternary systems 0, Q, 1^ , this choice of interactions does not 
apply. 

A second class of diagram occurs if the lipid cholesterol aligning interaction, J/c, is domi- 
nant. This produces a gel phase which is swollen with cholesterol. As there is no experimental 
evidence for this, we do not consider this choice of interactions further. 

A third class of diagram occurs when the interaction between cholesterols is at least 
comparable to that between lipids. Binary phase diagrams based on our calculations are 
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FIG. 2: Figure 2a. Binary phase diagrams for Ji^ = 0-78 Ju and Jcc = 0.73 Ju, and Ju as in Fig. 
1. The saturated Hpid-cholesterol mixture has a triple point very near the main-chain transition 
temperature, so that the gel, Id coexistence region is very narrow. Dashed lines are extrapolations. 
The ternary mixture at T = 300K is shown in Fig. 2b. 



shown in Fig. 2(a). In the cholesterol, saturated lipid system, there is coexistence of two 

liquid phases over a range of high temperatures, and a temperature of three-phase coexistence 

with a gel phase, consistent with refs A cut through the ternary diagram at a 

temperature of T = 300K below the triple temperature is shown in Fig. 2 (b). There 

is a large region of coexistence of the lo. Id, and gel phases. There are three regions of 

two-phase coexistence which extend to the binary axes. The coexistence region between lo 

and Id Uquids would be identified as the region in which rafts could exist. The diagram is 

□ 

topologically identical to that reported in ref j9|| . 

The degree of order, Sqd, in the lipid chains is given in terms of the second Legendre 
polynomial I^cdI = |-P2(^fc)|/2, where 6k is the angle between the normals to the bilayer and 
to the plane of the fc'th CH2 group. The order parameter for the saturated lipids is shown 
in Fig. 3 for each of the three phases which coexist in the phase diagram of Fig. 2b. As 
can be seen, the chains in the lo phase are indeed more ordered than those in the Id phase. 
This increased order is due to the interactions of the lipids with cholesterol, which is itself 
well ordered. We have also verified that the addition of cholesterol tends to disorder the gel 
phase, in agreement with experiment [l5|. 

The evolution of the ternary diagram with increasing temperature is simple. The region 
of three-phase coexistence shrinks and vanishes at the triple temperature of the binary 
cholesterol, hmp lipid system. Above the main chain temperature of this lipid, there is a two- 
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FIG. 3: The order parameters for the saturated tails in the three coexisting phases of the ternary 
mixture at 300K. Cholesterol concentration is greatest in the lo phase, which is more ordered 
than the Id phase. The effect of the larger cholesterol concentration is evident in the plateau 
corresponding to ordering about cholesterol's ring structure. 

phase region of lo and Id coexistence which stretches from one binary system with cholesterol 
to the other. At higher temperatures, this region detaches from the binary cholesterol. Imp 
lipid axis at a critical point, and the region shrinks with further increase of temperature 
until it vanishes at a critical point on the cholesterol, hmp lipid axis. It is clear from the 
evolution of our calculated ternary diagram that the liquid, liquid coexistence within our 
ternary diagram is intimately tied to a similar coexistence in the binary cholesterol, saturated 
lipid system. Given that the coexistence of liquids has been observed in the ternary systems, 
the calculation gives strong support to the view that such coexistence already exists in the 
cholesterol, hmp lipid system as well. 

We are grateful to Sarah Keller and Sarah Veatch for numerous profitable discussions. 
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0140500 and GTS-0338377. 
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